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ABSTRACT 

Array  weighting  designs  of  the  Dolph-Chebyshev  and  Kaiser-Bessel  type  are 
based  mathematically  on  orthogonal  polynomials.  The  theoretical  properties  of 
these  polynomials  give  rise  to  the  desirable  properties  of  the  resulting 
arrays.  This  paper  presents  results  for  array  weights  based  on  a  very 
general  set  of  orthogonal  polynomials  called  the  Jacobi  polynomials.  Many 
interesting  array  far-field  beampatterns  are  exhibited.  A  practical  means  of 
computing  all  the  array  weights  exactly  by  means  of  one  fast  Fourier  transform 
(FFT)  is  given.  This  method  is  quick  and  accurate  and  can  compute  the  weights 
for  arrays  having  large  numbers  of  elements.  It  can  efficiently  compute  both 
Dolph-Chebyshev  and  discrete  Kaiser-Bessel  weights  as  special  cases. 


ADMINISTRATIVE  INFORMATION 

This  technical  memorandum  was  prepared  under  NUSC  IR/IED  Project  No. 
A70212,  "Jacotri  Digital  Filter  and  Antenna  Design,"  Principal  Investigator,  R. 
L.  Streit  (Code  3232).  The  Program  Manager  is  CAPT.  Z.  L.  Newcomb,  Chief  of 
Naval  Material  (MAT  05B). 
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I.  INTRODUCTION 

Array  weights  of  the  Dolph-Chebyshev  and  Kai  ser-Bessel  type  are  based  on 
orthogonal  polynomials.  The  desirable  properties  of  these  arrays  are  due 
entirely  to  the  properties  of  the.underlyi  ng  orthogonal  polynomials.  In  this 
paper  the  mathematical  design  methodology  developed  in  [1],  which  parallels 
the  techniques  of  the  Dolph-Chebyshev  designs,  is  further  explored  for  Jacobi 
polyncmi  als. 

Analytical  expressions  for  the  underlying  weights  are  not  sought  here,  as 
they  are  in  [1],  because  such  expressi  ons  are  probably  not  competitive  with 
the  exact  FFT  based  method  presented  in  this  paper.  The  primary  purpose  of 
this  paper  is  to  present  a  unified  FFT  based  method  for  computing  array 
weights  based  on  any  of  the  Jacobi  polynomials.  This  method  is  quick  and 
accurate  and  can  compute  the  weights  for  arrays  having  large  numbers  of 
elements.  The  Appendix  gives  a  short  Fortran  program  for  computing  the 
weights  (given  a  subroutine  for  the  FFT).  This  program  can  efficiently 
compute  both  Dolph-Chebyshev  and  discrete  Kai  ser-Bessel  weights  as  special 
cases. 


This  paper  represents,  in  a  sense,  a  completion  of  certain  ideas  about 
array  design  using  orthogonal  polynomials.  Dolph  was  apparently  the  first  to 
use  an  orthogonal  polynomial  for  designing  an  array  weighting  function.  The 
polynomial  he  used  was  the  Chebyshev  polynomial  of  the  first  kind,  Tn(x), 
and  he  was  able  to  prove  an  optimality  ccndi  ti  on.  Unfortunately  his  proof  of 
this  optimality  ccndi  ti  on  relies  on  the  unique  behavi  or  of  the  graph  of 
Tn(x),  and  so  generalizati  ons  of  the  optimality  ccndi  ti  on  seem  unlikely. 
Nonetheless,  the  use  of  different  orthogonal  polynomials  can  lead*8seful 
weighting  functions.  For  example,  if  the  Chebyshev  polynomial  of  the  second 
kind,  Un(x),  is  used  in  place  of  Tq(x)  a  weighting  function  of 
Kai  ser-Bessel  type  results  (see  [1 J ) .  It  is  also  possible  to  use  a  general 

family  of  orthogonal  polynomials  that  contains  both  Tn(x)  and  Un(x)  as 
special  cases.  The  Gegenbauer  polynomials,  Cj((x),  are  one  such  family; 


that  is,  Tn(x)  and  Un(x)  are  the  special  cases  u  =  0  and  y  =  1, 
respectively.  This  family  is  used  in  [1]  and  gives  useful  and  interesting 
designs.  The  most  general  of  the  so-called  classical  orthogonal  polyncminals 
that  contains  all  these  examples  as  special  cases  is  the  family  of  Jacobi 
polynomials,  Pn  0  (x).  For  a  =  e  =  y  -  1/2  they  reduce  to  C^(x). 

Do  Jacobi  polynomials  turn  out  to  be  useful?  The  examples  presented  in  this 
paper  indicate  that,  although  many  interesting  new  array  designs  are  possible 
using  the  Jacobi  polyncminals,  the  most  useful  designs  in  this  general  family 
are  probably  those  that  have  already  been  discussed.  Consequently  the  new 


Jacobi  desi  gns 
problem. " 


might  be  said  to  be,  at  present,  "a  soluti  on  locking  for  a 
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II.  WEIGHT  GENERATION  BY  FFT 

The  far-field  beampattern  of  a  general  linear  beamformer  for  a  linear 
equispaced  array  having  2N  +  1  elements  (N  1)  with  element  positions  X|<  = 
k  x/(vD),  k  =  0,+  1,  ...,+  N,  is  given  by: 


N 

F(u)  =  Y1  w.  exp(-i2ir  ku/(vD))  (1) 

k  =  -  N  K 

where  the  integer  D  _>  1  is  given  and  v  >  0  is  a  fixed  real  constant  (vD  is  the 
number  of  elements  per  wavelength),  x  is  the  wavelength  of  the  design 
frequency,  and  u  is  defined  by 


u  =  sin  ea  -  sin  es. 


(2) 


where  ©s ,  — ir/2  &s  <it/2,  is  the  steering  (look)  angle  and  ea,  -  ir/2 
ea  <_ir/2,  is  the  arrival  angle  of  a  plane  wave.  Both  angles  are  measured 
from  a  line  normal  to  the  array  axis.  The  weights  (w^J  can  be,  in  general, 
any  set  of  complex  constants. 

Define  the  functions 


tDU)  - 

« 

pnu)  - 

n 

ZT 

Is 

,  n  >  0, 


(3) 

(4) 


where  ja|<jand  jb|<|  are  specified  constants.  By  simple  algebra 


Pn^D^z^ 

where  \c^\  depend  on  both 
gives 


nD 

k  =  -nD 


ckz 


(5) 


ja^l  andjb|<j  .  Substituting  z  =  exp(-iiru/(vD) ) 


(tD(exp(-i uu/  (vD) ) ) ) 

(6) 

nD 

, 

T]  c.  exp{-i2irku/(2vD) ). 

k  =  -nD  K 

(7) 
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By  comparing  (7)  with  (1),  it  is  seen  that  H(u)  is  the  far -field  beampattern 
of  a  linear  array  with  2nD  +  1  elements,  equispaced  x/(2vD)  apart,  and  with 
the  weight  applied  to  the  k-th  element.  It  will  be  shown  that  the  array 
weights  C|<  can  be  computed  from  function  values  of  tQ  and  Pn  by  means  of  an 

FFT. 


When  some  of  the  weights  C|<  =  0,  the  corresponding  elements  may  be 
eliminated  from  the  physical  array  without  altering  the  array's  far -field 
beampattern.  Elimination  of  zero  weighted  elements  (whenever  possible) 
minimizes  the  total  number  of  elements  required.  This  consideration  is 
especially  impoatant  when  tp  and  Pn  are  chosen  so  that  Pn(to(z))  is 
either  even  or  odd  in  z,  for  then  about  half  the  elements  need  not  be 
physically  present.  This  is  discussed  further  in  section  III. 

We  now  show  that  the  weights  jc^j.  in  (7)  can  be  computed  exactly  from 
(3)  and  (4)  by  the  fast  Fourier  transform  (FFT).  It  is  stressed  that  this 
procedure  is  theoretically  exact ,  not  approximate,  for  the  |c|<}  . 

Let  f(u)  be  any  complex  valued  function  of  a  real  variable  u  which  can  be 
written  exactly  in  the  form 


f(u)  =  £ 

k  -  -P 


d,  e 
k 


-iku 


(8) 


for  some  complex  constants  jd^}  .  Let 


2p  -  1 

f  _  1  2  F  e  2  trk  j  /  ( 2  p )  |<  _  q  1 

J  .  0  J  •  ’  . 


2P  -  1,  (9) 


where 


Fj  =  (-l)J  f((p-j)ir/p)  ,  j  =  0,  1 . 2p  -  1. 


(10) 


Thus  {f^l  is  the  inverse  FFT  of  order  2p  of  the  sequence  fFj}.  Substitute 
(8)  into  (10),  and  then  into  (9)  to  get 


2p  -  1 

i  E 

j  =  o 


l 

2P 


(-l)j  £ 

q  =  -p 


d  e 

q 


-i  irq  (p-j  )/p 


iirkj/p 


p 

q?-p 


(-l)q  d/  rt1  (-1)J  eMq  +k,J,p 


j  =  0 


k=0,l, . . . ,2p-l. 
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The  inner  sum  equals  2p  when  q  +  k  =  +  p,  +  3p,  ...  and  equals  zero 
otherwise.  Since  dq  =  0  for  jq)  >  p. 


-•f 


=  f (-l)p  (d  n  +  d  ),  if  k  =  0 


-P  P 


V, 


(-l)p_k  d 


p-k 


»  k— 1,  . ..,  2p— 1 . 


(11 -a ) 


(11. b) 


Now  let  f(u)  =  H(vDu/ir),  where  H(u)  is  given  by  (7),  and  p  =  nD.  Then 


Fj  =  (-l)J  Pn(tD(e-MnD  -  j)/nD)),  j  =  0,  1,  ...,  2nD-l, 


and  the  coefficients  C|<  in  (7)  are  just 


bn 


ck  "  « 


(-1) 


k-nD 


fnD-k 


,  if  k  =  -nD, 

,  if  k  =  -nD+1, 


,  if  k  =  nD, 


nD-1 


(12. a) 
(12. b) 

(12. c) 


where |fk^nD-1  is  the  inverse  FFT  l116  coefficients 

c_nD  and  cnD  cannot  be  computed  directly  by  FFT  because  of  aliasing,  as 

indicated  in  (11. a);  however,  by  direct  appeal  to  the  defining  equations  (3)  - 
(6),  it  follows  that  c_nn  and  cnQ  are  as  given  by  (12. a)  and  (12. c), 
respectively.  In  fact,  (11. a)  can  be  used  as  a  check  on  numerical  accuracy  in 
the  computations. 


It  should  be  obvious  that  some  special  forms  of  tp  and  Pn  can  be  used 
advantageously  to  reduce  the  size  of  the  FFT  required  to  compute  (12).  In 
general,  however,  no  special  structure  exists  and  the  smallest  FFT  size  that 
can  be  used  has  order  2nD. 


In  cases  where  2nD  is  not  an  integer  power  of  two,  the  FFT  is  still 
applicable  by  zero  filling  in  tp  and  Pn.  That  is,  D  is  replaced  by  the 
smallest  power  of  two  which  exceeds  or  equals  D,  say  D1.  The  function  tg  is 
then  merely  considered  to  be  a  special  case  of  tg».  Similarly,  Pn  is  a 
special  case  of  Pn •  for  some  smallest  power  of  two,  n1,  which  is  greater 
than  or  equal  to  n.  The  required  size  of  the  FFT  is  thus  2n ' D 1 .  The 

coefficients  {  C|<}  are  still  given  by  (12);  however,  from  (7),  it  must  be 
the  case  that  C|<  =  0  forjk  |>  nD. 

Hie  Appendix  gives  a  Fortran  subroutine  for  computing  the  array  weights 
by  an  FFT,  given  subroutines  for  evaluating  Pn(z)  and  tg(z).  The  program 
is  specialized  to  the  case  D  =  1,  but  it  can  be  easily  altered  to  accommodate 
larger  values  of  D.  The  program  is  also  written  on  the  assumption  that  2nD  is 
a  power  of  two.  As  discussed  in  the  preceding  paragraph,  zero  filling  allows 
the  most  general  situation  to  go  through. 
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III.  THE  TEN  PARAMETER  JACOBI  WEIGHT  FAMILY 


The  most  important  special  case  of  H(u)  seems  to  be  D  =  1.  Although 
there  may  be  interesting  possibilities  when  D  >  1  (consider,  for  example,  the 
identity  Tn(TD(cos  e))  =  Tn+n(cos  e),  where  Tn  is  the  Chebyshev 
polynomial  of  the  first  kind),  these  cases  are  not  explored  in  this  paper. 


Before  proceeding,  it  is  instructive  to  see  first  how  the  usual 
Dolph-Chebyshev  case  for  half  wavelength  equispaced  arrays  is  derived  as  a 
special  case  of  H(u).  Let 


N  =  number  of  array  elements 


tj(z)  =  Zq(z-1  +  z)/2 
Pfi(z)  =  Tpf(z) , 


where  Tff(z)  is  the  Chebyshev  polynomial  of  the  first  kind  and  the  real 
constant  Zq  is  given  by 


Z0  "  H(Q  +  (q2_1)  1/2)1/"  +  (Q  -  (q2  -1)1/2)1/"J  (15) 


where 


Q  =  10  |  S  | / 20 

S  =  specified  sidelobe  level  (in  dB). 

For  these  values  it  follows  that  ti  (exp(-iiru/2))  =  Zn  cos(iru/2)  and,  from 
(6)  and  (7), 


H(u)  =  Tfi(  ZQ  cos(irU/c ) )  (16) 


n 

=  Z  ck  exp(-i wku/2 ) .  (17) 

k  =  -n  K 

By  comparing  (17)  to  (1),  it  would  appear  at  first  glance  that  this  array  is 
quarter  wavelength  equispaced  with  2n  +  1  =  2fi  -  1  elements.  However,  every 
other  coefficient  in  (17)  is  identically  zero  because  Tff(z)  is  always  either 
an  even  or  an  odd  function  in  z.  Deleting  the  zero  weighted  elements  reduces 
(17)  to  two  slightly  different  cases,  depending  only  on  whether  ft  is  even  or 
odd.  These  cases  are  not  given  explicitly  here.  Note,  however,  that  ft  even 
implies  that  ft  is  odd  and  that  Tfi’(z)  contains  no  even  powers  of  z,  which 
means  that  Tff(z)  has  (ff  +  1 ) / 2  nojn-zero  coefficients  and,  consequently,  that 
Tfi(ti(z))  has  precisely  fT  +  1  =  N  non-zero  terms  in  the  expansion  (17). 
Simitar  reasoning  holds  for  tt  odd. 


9 


TM  851015 


To  summarize:  The  Dolph-Chebyshev  array  design  is  thought  of,  in  the 
context  of  this  paper,  as  a  quarter  wavelength  equispaced  array  in  which  half 
the  elements  (every  other  one)  has  been  zero  weighted.  Dolph-Chebyshev  arrays 
of  both  even  and  odd  numbers  of  elements  are  thought  of  in  this  way. 

From  (3),  the  most  general  form  for  D  =  1  is 

t^(z)  =  a_i  z-1  +  ao  +  ai  z. 

For  reasons  that  will  become  clear,  the  slightly  more  restrictive  form 


tx (z )  =  ZqO""1  z'1  +  aQ  +  rQ  z)/2,  r  ±  0,  (18) 

is  adopted,  where  zo,  rQ,  and  ao  are  arbitrary  complex  constants.  The 
only  useful  form  excluded  by  (18)  is  obtained  by  changing  the  sign  of  the 
highest  order  term;  this  latter  form  corresponds  to  "difference"  patterns  and, 
for  ease  of  exposition,  is  not  discussed  further.  For  the  remainder  of  the 
paper,  (18)  is  taken  as  the  definition  of  ti(z).  Note  that  tife-i0-)  = 
zq  cos  ©•  if  ro  =  1  and  ag  =  0. 

The  most  general  class  of  polynomials  Pn(z)  considered  in  this  paper  are 
the  Jacobi  polynomials,  denoted  P^<*>S)(z).  They  are  defined  explicitly  by 


>nMlz)  'ST  k?0(k)(nW1)k  <“',ktl>n4c(Lii)k  <19> 

for  all  complex  values  of  a,  e,  and  z.  For  a  >  -1  and  a  >  -1,  the  Jacobi 
polynomials  are  orthogonal  on  the  real  interval  -1  <  z  <_ +  1,  but  they  are  not 
necessarily  orthogonal  for  other  values  of  a  and  b.  The  best  available  method 
for  computing  Pn(<*,B)(z)  relies  not  on  (19)  but  on  the  three  term 
recursion  which  they  satisfy.  The  published  algorithm  [2],  [3]  based  on  this 
recursion  is  easily  modified  to  compute  the  Jacobi  polynomials  for  complex  a, 
B,  and  z  for  all  values  of  the  degree  n  that  are  likely  to  be  of  practical 
interest,  say  n  <  150.  A  thorough  mathematical  treatment  of  Pn(“>6)(z) 
is  available  in  [4]. 

The  generalized  Laguerre  and  Hermite  polynomials,  together  with  the 
Jacobi  polynomials,  constitute  a  complete  list  of  the  so-called  classical 
orthogonal  polynomials.  Array  designs  can  also  be  based  on  the  Laguerre  and 
Hermite  polynomials  and  will,  of  course,  be  different  from  those  based  on 
Jacobi's.  Although  these  designs  are  probably  interesting,  in  this  paper 
attention  is  restricted  to  the  Jacobi  polynomials. 

The  use  of  (18)  and  (19)  givesrise  to  a  five  parameter  family  of  weights 
which  includes  nearly  all  the  well  known  analytic  families  of  weights  as 
special  cases.  (The  most  prominent  exception  is  Taylor  weighting).  The  five 
parameters  are  zq,  ro,  ag,  a,  and  b>  Each  parameter  can  be  complex,  so 
there  are  actually  10  real  parameters  if  the  real  and  imaginary  parts  of  each 
are  counted  separately.  The  Fortran  program  listed  in  the  Appendix  is  written 
for  this  general  case. 
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The  simplest  way  to  explore  the  properties  of  these  ten  parameters  is  to 
perturb  each  parameter  separately  while  holding  the  others  fixed  at  some  nominal 
value.  The  nominal  parameter  values  chosen  here  for  the  examples  are  those  that 
give  rise  to  the  Dolph-Chebyshev  design  for  an  array  of  33  elements  with  a 
sidelobe  level  of  -30  dB.  Specifically, 

aO  =  -.50  +  0.0  i 
SO  =  -.50  +  0.0  i 
ag  =  0.0  +  0.0  i 

ro  =  1.0  +  0.0  i  (20) 

ZQ  =  Zo  =  1.008408  +  0.0  i, 

where  Zq  is  computed  from  (15)  with  S  =  -30  dB  and  n  =  32.  Recalling  earlier 
remarks  in  this  section  concerning  the  interpretation  of  half  wavelength 
equispaced  arrays  as  quarterwavelength  equispaced  arrays  with  zero  weights,  set 
n  =  32  in  (6)  and  (7).  Thus,  in  principle,  the  example  is  a  2n  +  1  =  65  element 
quarterwavelength  equispaced  array.  The  coefficients  Ecu]  are  computed  from 
(12).  In  (12. a)  ana  (12. c),  note  that  the  identity  [4,  Eq.  (4.21.6)] 


holds  and  so,  from  (18), 


a0  -  f  z0r0’  a-D  -  ?  z0r0_1  <22> 

Each  of  the  ten  parameters  is  both  increased  as  well  as  decreased  from 
its  nominal  value  given  in  (20).  Thus  there  are  21  cases,  including  the 
nominal  case  itself  in  (20).  table  1  displays  these  cases  and  gives  each  an 
identifying  case  name.  To  each  case  in  Table  1  there  corresponds  a  graph  of 
the  far-field  beampattern  H(u)  on  the  u-interval  [0,4]  and  two  bar  charts,  one 
for  the  real  part  and  one  for  the  imaginary  part  of  the  weights  corresponding 
to  that  case. 


Figure 

Number  Case  Names  Value  of  Perturbed  Parameter 


1  NOW  no  deviations  from  (20) 


2 

Z.l 

Z.2 

Z.3 

Z.4 

Z0+.003 

Z0-.003 

zo+- 003i 

zo-.003i 

3 

A.l 

A. 2 

A. 3 

A. 4 

ao+.003 

ao-.003 

ao+. 003i 

ao-.003i 

4 

R.l 

R.2 

R.3 

R.4 

r0+.03 

ro-«  03 

ro+.03i 

ro-.03i 

5 

a.l 

a. 2 

a. 3 

a. 4 

<*0+*3 

a0-«3 

a0+*3i 

a0-.3i 

6 

e.l 

B.2 

e.3 

B  .4 

Bq+-3 

BQ-.3 

BQ+.3i 

BQ-.3i 

Table  1  Perturbed  parameter  values;  deviation  from  the  nominal  values  (20) 

In  general,  H(u)  is  periodic  with  a  period  of  length  2vD.  Since  D  =  1 
and  v  =  2  in  these  examples,  any  interval  of  length  4  suffices  to  exhibit  all 
the  structure  of  H(u). 
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In  all  the  bar  charts  presented,  upward  lines  indicate  positive  weights 
and  length  is  proportional  to  magnitude.  Similarly,  downward  lines  indicate 
negative  weights.  These  upward  and  downward  lines  are  ordered  from  left  to 
right  and  correspond  to  elements  numbered  from  -32  to  +32.  Any  element 
receiving  a  zero  weight  is  indicated  by  a  simple  "x"  marking  its  position.  In 
particular,  notice  that  the  nominal  case,  NOM,  of  Figure  1  has  only  33 
non-zero  weights.  The  nominal  case,  as  has  been  said,  is  a  half  wavelength 
equispaced  array  being  treated  as  a  quarter  wavelength  equispaced  array.  The 
weights  in  each  case  are  normalized  by  the  largest  magnitude  of  any  real  or 
imaginary  part;  thus,  the  normalization  between  cases  is  not  exactly  the  same. 

Figure  1  is  the  reference  case  (20)  and  needs  no  further  comment. 

Figure  2  perturbs  only  z q.  The  cases  Z.l  and  Z.2  are  expected  since 
zq  is  merely  increased  or  decreased  in  its  real  part  alone.  When  zo  is 
perturbed  by  adding  an  imaginary  component,  the  array  still  has  33  non-zero 
weights  and  so  is,  in  effect,  half  wavelength  equispaced.  It  is  surprising 
how  much  can  be  added  to  the  imaginary  parts  of  the  weights  without  seriously 
degrading  the  beampattern.  The  beampatterns  in  Z.3  and  Z.4  are  identical. 

Figure  3  perturbs  aQ  from  its  nominal  zero  value.  Any  perturbation 
produces  a  quarter  wavelength  equispaced  array  which  is  symmetrically 
weighted.  One  way  to  discuss  the  results  is  to  visualize  the  65  element  array 
as  being  composed  of  two  half  wavelength  equispaced  arrays — one  having  33 
elements  and  the  other  32  elements  with  the  elements  of  the  two  arrays 
interlaced.  Thus,  perturbing  the  real  part  of  ao  is  equivalent  to  adding  or 
subtracting  the  outputs  of  these  two  arrays.  Perturbing  the  imaginary  part  of 
ag  is  equivalent  to  adding  or  subtracting  the  outputs  after  first  putting 
them  in  phase  quadrature  with  respect  to  each  other.  The  beampatterns  A. 3  and 
A. 4  are  identical  to  each  other,  but  they  are  NOT  the  same  as  Z.3  and  Z.4. 

Figure  4  perturbs  ro  from  its  nominal  value  of  +1.  Any  small 
perturbation  produces  asymmetrically  weighted  half  wavelength  equispaced 
arrays.  Real  perturbations  of  ro  produce  only  real  weights  and  have 
beampatterns  without  any  true  nulls.  Pure  imaginary  perturbations  do  not 
alter  the  beampattern  from  its  nominal  case,  even  though  the  weights  develop 
an  interesting  sinusoidal  character  in  their  imaginary  parts. 

Figure  5  perturbs  a  from  its  nominal  value  of  «o  =  -1/2.  Any 
perturbation  produces  a  quarter  wavelength  equispaced  array  which  is 
symmetrically  weighted.  Real  perturbations  yield  real  weights  while  pure 
imaginary  perturbations  yield  complex  weights.  The  first  grating  lobe  in  case 
o.l  is  at  about  -8  dB  instead  of  0  dB;  the  same  is  true  of  the  MRA  in  case  a. 2. 

Figure  6  perturbs  b  from  its  nominal  value  of  bo  =  -1/2.  Any 

perturbation  produces  quarter  wavelength  equispaced  arrays  which  are 
symmetrically  weighted.  Real  perturbations  yield  real  weights,  while  pure 
imaginary  perturbations  yield  complex  weights.  The  first  grating  lobe  in  case 
3.2  is  suppressed  to  about  -8  dB,  while  in  case  b.1  the  MRA  is  depressed  to 

-8  dB.  Figure  6  and  Figure  5  should  be  compared. 
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The  following  observations  seem  to  hold: 

1.  The  real  part  of  a  controls  the  "upper"  envelope  of  H(u)  near  the 
center  peak . 

2.  The  absolute  value  of  the  imaginary  part  of  a  controls  the  "lower" 
envelope  of  H(u)  near  the  center  peak . 

3.  The  real  part  of  e  controls  the  "upper"  envelope  of  H(u)  near  the 
first  grating  lobe . 

4.  The  absolute  value  of  the  imaginary  part  of  e  controls  the  "lower" 
envelope  of  H(u)  near  the  first  grating  lobe. 

Other  parameters  (the  imaginary  parts  of  a  and  zq)  also  affect  the  "lower" 
envelope  of  H(u),  but  the  dominant  effects  seem  to  be  due  to  the  imaginary 
parts  of  o  and  $•  The  imaginary  part  of  r  does  not  affect  the  lower  envelope 
at  al 1 . 

By  changing  the  parameters  simultaneously  in  different  ways,  the 
different  effects  may  be  combined,  at  least  for  small  perturbations.  Examples 
of  this  are  not  included  here. 
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Figure  1.  No  perturbations;  the  nominal  case  (20) 
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Figure  2.  Perturbations  of  z 
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IV.  SUMMARY  AND  CONCLUDING  REMARKS 

It  has  been  shown  that  array  weights  based  on  the  Jacobi  orthogonal 
polynomials  can  be  computed  exactly  by  means  of  FFT.  As  a  special  case, 
weights  based  on  the  Gegenbauer  polynomials  can  also  be  computed  exactly  by 
FFT,  instead  of  analytically  as  in  [1].  Examples  have  been  presented  to  show 
the  effects  of  varying  the  ten  parameters  in  the  Jacobi  family. 

Further  work  in  this  area  is  possible.  In  addition  to  the  Jacobi 
polynomials,  one  may  also  use  the  generalized  Laguerre  and  the  Hermite 
polynomials.  In  fact,  any  orthogonal  polynomial  family  that  has  interesting 
structural  features  can  be  the  basis  of  a  weighting  family  which  inherits  this 
structure.  In  a  different  direction,  certain  cases  for  D  >  1  may  yield 
interesting  designs  and  have  not  been  explored.  The  weights  corresponding  to 
all  these  cases  can  computed  exactly  by  the  FFT  method  presented  in  this  paper. 
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APPENDIX 

The  Fortran  program  JACWTS  listed  below  assumes  that  D  =  1  and  that  2nD 
is  a  power  of  2.  The  function  to(z)  is  defined  exactly  as  in  (18),  and  the 
polynomial  Pn(z)  is  taken  to  be  the  Jacobi  polynomial  Pn(<*>s)(z). 

This  program  is  an  implementation  of  the  (exact)  FFT  method  described  by 
Eq.  (12),  where  a±g  and  bn  are  given  by  (22)  and  (21),  respectively.  The 
user  of  JACWTS  need  only  specify  values  for  a,  B,  ag,  rg,  and  z q.  In 
JACWTS  these  variables  are  referred  to  by  the  labels  ALPHA,  BETA,  AO,  RO, 

ZO,  respectively.  The  arrays  X  and  Y  contain,  on  output,  the  real  and 
imaginary  parts  of  the  array  weights  [c^}  .  These  two  arrays  must  be 
dimensioned  at  least  2nD+l  in  the  routine  which  calls  JACWTS.  The  integers  n 
and  D  are  referred  to  by  the  labels  N  and  D,  respectively,  in  the  subroutine 
argument  list.  Also,  LOGN  and  LOGO  are  defined  so  that  N  =  2**LQGN  and  D  = 
2**L0GD. 

This  program  assumes  that  a  subroutine  named  JACOBI  evaluates  the  Jacobi 
polynomial  (19)  for  arbitrary  complex  values  of  a,  s,  and  z.  This  subroutine 
can  be  based  on  the  published  codes  in  [2]  and  [3].  This  program  also  assumes 
that  subroutines  are  available  for  computing  a  complex  FFT  of  size  2nD;  the 
particular  ones  used  here  are  based  on  Markel 's  method  and  are  not  listed. 
Their  names  are  DPMCOS  and  DPMFFT.  These  routines  require  a  work  array,  C, 
dimensioned  at  least  2nD  in  the  routine  which  calls  JACWTS. 

JACWTS  is  written  in  double  precision  complex  mode  to  forestall  any 
numerical  round-off  error  problems  that  might  arise.  The  test  suggested  in 
Section  II  (that  follows  from  the  resolution  of  the  aliasing  effects  as  in 
(12))  is  incorporated.  It  is  the  only  test  used  to  ascertain  whether 
numerical  round-off  of  significant  proportions  occurred.  No  numerical 
difficulties  have  been  detected  by  this  test  to  date,  which  indicates  that  the 
computation  is  usually  numerically  reliable. 
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SUBROUTINE  JACWTS  (  X ,  Y ,  C ,  N ,  LOGN  ,  D ,  LOGO ,  ALPHA ,  BETA  ,  AO  f  RO  ,  ZO  ) 
COMPLEX* 16  Z, T, H,  S, R ,TSUBO,  ALPHA, BETA, ZO, AO, RO, JACOBI 
INTEGER  M, LOGN,D, LOGO, TWO ND 

DOUBLE  PRECISION  EPSI , ARG , P I , X C 1) , Y ( 1 ) , C C 1 ) 

OATA  PI, EP3I/3. 14159265358979 3 23800,0.50-7/ 

M=LQGN  +  LOGD  +  1 

ND=N*D 

TViONO*2*HO 

...  ALL  WEIGHTS  EXCEPT  THE  FIRST  AMD  THE  LAST 
I=+l 

00  10  J=t , TWQND 
AflG  =  -PI*(.MD-J+U/.MD 
Z-OCMPLX(COSCARG) , SIN (ARG)) 

T=TSUSO(0,AO,RO,ZO,Z) 

HsJACOBKNf  ALPHA,  3EIA,T) 

X(J)=I*QR2AL(H) 

YCJ)  =  I*DUAG(H) 

I=-I 

CONTINUE 

CALL  DPMCQS  (C,  T'WOND) 

CALL  DPMFFT(X,Y,C,Mt+l) 

Tati 

00  15  J=1,TNQND 
X(  J)  =I*X(J) /TWQrJQ 
YCJ)aI*Y(J)/TI«OND 
t  =  -I 

CONTINUE 

...  THE  FIRST  AND  LAST  WEIGHTS 

H=.25D0*Z0*R0 

S=.25D0*Z0/R0 

R  =  1 « ODO 

T  =  l ,000 

DO  18  1=1, M 

Z=(  (M+il-Itl  )+ALPHA  +  3ETA)/I 

T=T*H*Z 

R=R*S*Z 

CONTINUE 

X(1)=0RSALCR) 

Y(l)=OIMAG(R) 

X(TWQND+n=OREALCT) 

Y(TW0ND+1)=DIMAGCT) 

...  NUMERICAL  ACCURACY  TEST 
ARG=(A3S(X(t  )-OR£AL(R  +  T)  )  tABSCYCn-DIMAGCR+T)  )  ) 

+  /(l.ODO  +  AftSC <(1))  +  A3S(Y(1)  )  ) 

IFI ARG.GT.EPSDPRIMT  50 

FORM Alt  *  NUMERICAL  ROUND-OFF  ERROR  IS  SIGNIFICANT. ') 

RETURN 

END 

FUNCTION  TSUBDCD, AO,RO,ZO,Z) 

COMPLEX* 1 6  Z , ZO , AO , RO , TSU80 
INTEGER  D 

TREAT  THE  CASE  0=1?  IGNORE  OTHER  VALUES. 

TSU3D=.500*Z0*(  (1 . OPO / ( RQ*Z ) )  +  AO  +  (RO*Z)  ) 

RETURN 

END 
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